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Abstract. - We investigate the inherent structure (IS) dynamics of mean-field fimte-size 
spin-glass models whose high-temperature dynamics is described in the thermodynamic limit 
by the schematic Mode Coupling Theory for super-cooled liquids. Near the threshold energy the 
dynamics is ruled by activated processes which induce a logarithmic slow relaxation. We show 
the presence of aging in both the IS correlation and integrated response functions and check 
the validity of the one-step replica symmetry breaking scenario in the presence of activated 
processes. Our work shows; 1) The violation of the fluctuation-dissipation theorem is given 
by the configurational entropy, 2) The intermediate time regime (log(t) ~ A'^) in mean-field 
theory automatically includes activated processes opening the way to analytically investigate 
activated processes by computing corrections beyond mean-field. 



After many years of research on the structural glass problem a lot of experimental data has 
been collected but a convincing theory is still needed jl], || . The two most successful theories 
for the glass transition are the Adam-Gibbs-DiMarzio and the ideal Mode Coupling Theory 
(MCT) Despite their different character both are mean-field theories. In the former case 
the mean-field aspect lies in the notion of configurational entropy which assumes a breaking 
of the phase space in disconnected ergodic components. In the latter the presence of the MCT 
transition marks the onset of ergodicity breaking where the diverging of a characteristic time 
occurs. Both approaches have been successfully unified in the context of spin-glass theories Q. 

To go beyond mean-field it is necessary to include activated processes, a very difficult task 
since it implies the knowledge of the excitations involved in the dynamics. Recent theoretical 
and numerical results clearly show that the slowing down of the dynamics near the structural 
glass transition is strongly connected to the complex topology of the potential energy landscape 
[^. In a glass- forming systems this is made by many deep valleys connected by saddles. The 
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time evolution can then be divided into an intra-valley and an inter-valley motion. When 
the temperature is lowered down to the order of the critical MCT temperature Tmct the two 
motions become well separated in time and relaxation dynamics, dominated by inter-valley 
processes, slows down displaying non-exponential behavior. 

To deal with this picture Stillinger and Weber (SW) j|] introduced the concept of inherent 
structure (IS) defined as the local stable minima of the potential energy reached through a 
steepest descent energy minimization process. All configurations which under this mapping 
flow to the same IS define the basin (of attraction) of the IS. It is now a simple matter to 
define an IS-based thermodynamics by replacing the partition sum with a sum over IS ^,|^. 
Direct consequence of this is the introduction of a configurational entropy Sc{e) which counts 
the number of different IS with the same energy e: Q.{e) — exp(A^Sc(e)). In is clear that as 
long as the configurations retained in an IS-based partition sum are those dominating the 
thermodynamics at the given temperature T, the equilibrium behavior is correctly described. 
However, despite of the fact that the Sc(e) so defined is a dynamical quantity, it is far from 
obvious that it describes the long-time non-equilibrium behavior. The reduction from the 
real dynamics to an IS dynamics is what, in the theory of dynamical systems, is called a 
symbolic dynamics. This describes correctly the dynamics only if it is associated to a generating 
partition |^ . In general for a generic dynamics it is not at all trivial to demonstrate that such 
a partition exists, and even if it does exist, how to find it. Nevertheless, we can argue that in a 
many valley dynamics with activated dynamics the SW mapping should be a "good" mapping. 
Indeed, since the SW mapping replaces each configuration in a IS-basin with the IS itself, it 
is clear it will be a good mapping for the long time non-equilibrium dynamics if the system 
spends a lot of time inside the basins. Under this assumption the dynamics on time scales 
larger than the typical residence time inside a IS-basin should be quite well described by the IS 
dynamics. Recent numerical results on Lennard-Jones mixtures |^ supports this scenario for 
supercooled liquids. However, this mapping may not be valid for all kind of non-equilibirum 
systems [|lO|. Actually, for coarsening systems the SW mapping is not meaningful | pl] |. 

In this Letter we extend the analysis to finite-size mean-field glass models and propose that 
activated processes seen in supercooled liquids can be treated at a mean-field level by including 
finite-size effects in the dynamics of an infinite mean-field system going beyond the saddle- 
point approximation, i.e., beyond the ideal MCT. This observation is quite reminiscent of the 
dynamical approach of Sompolinsky and opens the way to address activated processes 
in structural glasses starting from mean-field theories. After the work by Sompolinsky the 
inclusion of activated processes from this point of view has not attracted much attention [p^ 
and mean-field dynamical studies have mainly considered the N ^ oo limit before the large- 
time limit Other possible approaches analyze activated process by considering instanton 
solutions of the mean-field equations [p^ , [l5|] . 

This work is the natural continuation of a previous one where we introduced finite- 
size mean-field glasses for the analysis of the SW configurational entropy. Here we extend 
our study to the long-time non-equilibrium dynamics finding that, once finite-size effects are 
included and activated processes appear in a natural way, relaxational dynamics is driven by 
the configurational entropy giving support to the one-step replica-symmetry breaking (RSB) 
scenario beyond mean-field. The present results are in agreement with recent simulations on 
Lennard-Jones systems [ pT[ although here we go further and verify the explicit connection 
between the violation of the fluctuation-dissipation theorem and the conflgurational entropy. 
Moreover, our results give further support to the relevance of p-spin-like models for the de- 
scription of the glass transition in structural glasses. 

Following Ref. we consider the dynamics of the Ising-spin Random Orthogonal Model 
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(ROM) KAhm defined by the Hamiltonian 



= -2^J,,a,a, (1) 



where at = ±1 are N Ising spin variables, and Jij is a N x N symmetric random orthogonal 
matrix with Ju = 0. We use the heat-bath Monte Carlo scheme with random sequential 
spin updating. This model presents a thermodynamic glass transition typical of mean field 
p-spin glasses with p > 2 ||2^ at Tmct = 0.536 with threshold energy per spin eth = —1-87. 
Below this temperature the system is dynamically confined into a metastable state (basin 
or valley) and cannot reach true equilibrium. A equilibrium transition with collapse of the 
configurational entropy takes place below Tmct at the Kauzmann temperature Tc = 0.25, 
with critical energy per spin Cc = —1.936 The analysis of the free-energy landscape 



of this model |19| reveals that the phase space is composed by an exponentially large (in TV) 
number of different basins, separated by infinitely large (for N — > oo) barriers. Above Tmct 
and in the large iV limit the IS with e = eth attract most (exponentially in N) of the states and 
dominate the behavior of the system. For finite N basins of IS with e 7^ eth have statistical 
weight and may influence the dynamics |l6|| . 

To study the non-equilibrium relaxational dynamics we quench the system at time t — 
from an equilibrium state at temperature Tj > Tg to a final temperature Tf < Tg. The glass 
transition temperature Tg is defined, in accordance with the "experimental" definition, as the 
temperature below which we cannot equilibrate the system on the longest Monte Carlo run. 
The associated IS dynamics is obtained by regularly quenching the system down to T = 
from the relaxing configuration and recording the IS associated with the instantaneous basin. 
We shall consider both one-time quantities, as the average IS energy, and two-time quantities, 
as correlation and response functions. 

Relaxation of one-time quantities. In figure |l| it is shown the average IS energy per spin 
(e)(i) as function of time for a system with N — 300 spins. The analysis of the figure reveals 
that the relaxation process can be divided into two different regimes where {e){t) decreases 
with different power laws. The slope of the decay in the first regime is independent of Tf while 
the slope in the second regime is independent of both Tj and Tf [cfr. panel (a)]. For a given 
Tj the final temperature Tf only controls the cross-over between the two regimes. A similar 
behavior has been observed in molecular dynamics simulations of supercooled liquids @ . We 
note that since we use discrete variables, and hence a faster dynamics, the very-early regime 
observed in where (e}(t) is almost independent of t is absent. Moreover since we use fully 
connected systems the power law exponents change with N, even if the qualitative scenario 
is unchanged. The two regimes are associated with different relaxation processes. In the first 
part the system has enough energy and relaxation is mainly due to path search out of basins 
through saddles of energy lower than fceTf. During this process the system explores deeper 
and deeper valleys while decreasing its energy. The process stops when all barrier heights 
become of ©(fceTf ). From now on the relaxation can only proceed via activated processes and 
the dynamics slows down becoming logarithmic in time. Note that, for finite N, the activated 
regime starts already above the threshold indicating that not only stable states exist above 
eth but these also influence the hopping dynamics. 

Aging in correlation and response functions. More information on the non-equilibrium 
activated regime can be obtained from the analysis of correlation and response functions. In 
order to study the relevance of the SW mapping we consider the IS-based correlation and 
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Fig. 1 - Average IS as a function of time, (a): Filled symbols Ti = 3.0, empty symbols Ti = 0.9. 
Final temperatures are as in panel (6). (6): Ti — 3.0. The average is over different equilibrium initial 
configurations. 



response functions defined as, 

1 ^ 

Cis(t,s) = -^a,"(i)a°(.) (2) 

1=1 

1 ^ Sa^it) 

i—l ^ ^ 

where cr*'(t) is the IS at time t and hi an external field. It can be shown that in equilibrium 
the fluctuation-dissipation theorem TRis{t — s) = —dtCis{t — s) holds [^. To study the 
non-equilibrium response function we quench the system at time t = from an equilibrium 
state at Ti > Tg to T{ < Tg and, after a given waiting time t^, we make a clone of it to which a 
(small) constant uniform magnetic field is applied. To reduce fluctuations the same stochastic 
noise is used for both replicas. The experiment is repeated for different (small) fields to control 
the linear response regime. Results for the correlation function and the integrated response, 
or zero-field cooled susceptibility, xzFc{t,tw) — /j dt' Ris{t,t') are reported in figure ||. The 
presence of aging is rather clear. 

The one-step RSB scenario and the SW configurational entropy. A way to see how equipar- 
tition is broken in the non-equilibrium regime is through the fluctuation-dissipation ratio, 
initially introduced in the context of spherical models for spin glasses H, ^|2^, which 
gives a direct measure of frozen degrees of freedom in the response of the system. Usu- 
ally the broken of equipartition is expressed through an effective temperature deflned as 
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Fig. 2 - IS correlation (above) and integrated response (below) functions as function of time for 
different waiting times. The system size is A*' = 300, Ti — 3, Tf = 0.2 and Tg ~ 0.5. Data have been 
averaged over about 400 dynamical histories. 



Teff(i, iw) = dt„Cis{t,tvj)/Ris{t,tvj), or alternatively as 

rr-irn \ dxzFc{t,t-w) 

T..iC,s)-- g^^^^^^^^^ ■ (4) 

which is the slope of the curves xzfc versus Cis . This temperature reduces to T in equilibrium 
and is larger when equipartition is broken. 

In the context of mean-field theories for the glass transition T^ff is given by the constant 
temperature derivative of the free-energy with respect to the configurational entropy: = 
[dsc{f)/df]T @,|2^, which in the present IS analysis means 
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(5) 



where /is is the free energy of the IS with energy e, i.e., obtained from the partition sum 
restricted to the configurations in the IS-basins of IS with energy e. Writing fis — e + 
Ae(T, e) — Tsis{T,e) we see that [9/is(e)/<9e]T = 1 -I- i5(e,T), where the last term measure 
how IS with different e differ one each other. 

Results from our numerical simulations are shown in figure ^. For each t„ two distinct 
behaviors are clearly seen. For large values of Cis the slope of all the curves is /3f — 1/Tf. 
Below a fw-dependent value of Cjs the curve changes slope and the effective temperature 
increases. In order to compare our results with prediction (^, we have to compute /is. 
This can be done from the knowledge of the probability that an equilibrium configuration 
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Fig. 3 - Integrated response function as a function of IS correlation function. The data are from Fig. 
|2| The dash line has slope /3f — 5.0, while the full lines is the prediction ^ with 5{e, T) — and 
Sc(e) from Ref. reff(2") ~ 0.694, Tcft(2^^) ~ 0.634 and T^ft{2^^) ~ 0.608. The dot-dashed line 

is /Jeff for t„ = 2 drawn for comparison. 



at temperature T — \/ (3 lies in a basin associated with IS of energy between e and e -I- de: 
PN{e,T) = exp N\-I3fis + Sc(e) + Pf{T)] where f{T) is the full free-energy. Using the 
results from Refs. |l^,|l^] both /is and its derivative can be estimated. In all cases we found 
S{e,T) ~ 0, within numerical errors, a results in agreement with the fact that the ROM is 
rather similar to the Random Energy Model |^ for which S{e,T) — 0. Using for Sc(e) the 
expression obtained in Ref. ||l^ we obtain from (||) the slopes PcS = l/7cff shown in figure ^. 
The agreement is rather good. Data for Tf = 0.1 and 0.3 are more noisy, but consistent with 
this identification. 

We therefore conclude that in the activated regime the elfective temperature T^g derived 
from the SW configurational entropy agree extremely well with the numerical data confirming 
the one-step scenario. This is a highly non-trivial result. First of all, as discussed above, it is 
not obvious that the SW mapping correctly describes the long-time non-equilibrium dynamics. 
It is plausible that it works for an activated dynamics, but plausibility is not a proof. Second, 
equation is derived within mean-field, i.e., for TV —^ oo. Here we apply it for states which 
exist only for finite N being, as discussed above, those which govern the relaxational dynamics 
in this regime but disappear in the thermodynamic limit. Thus our finding gives to (|^) a 
broader validity. We note that quite new results from numerical simulations of Lennard- Jones 
mixtures are in agreement with (||), but with S{e,T) ^ p^ . 

From the discussion on the SW approach it follows that if the SW mapping defines a good 
mapping then it should equally well describe the long-time relaxational dynamics regardless 
which configuration we use to identify the IS-basin. To check this point we have repeated the 
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simulations using slightly different definitions of IS, i.e., changing the minimization rule. In 
all cases we have found the same results, so that in this case the SW mapping provides a good 
"symbolic" dynamics for the long-time non-equilibrium behavior. 

We now summarize our findings. The Stillinger- Weber decomposition of phase space in 
terms of IS is a natural and simple statistical description of a dynamical system with activated 
dynamics among different basins when the time scales for motions inside a basin and between 
basins are well separated. Moreover we have shown that finite-size mean-field glasses are 
valuable models to describe activated processes in structural glasses. In the activated regime 
these models display slow logarithmic relaxation with two class of motions which confirm the 
one-step replica symmetry breaking scenario: a fast intra-basin motion where fluctuation- 
dissipation holds and a slow inter-basin motion corresponding to activated processes. The 
fluctuation-dissipation ratio in the activated regime can be fully described in terms of the 
configurational entropy generalizing the mean-field scenario with the inclusion of activated 
processes. Our results suggest the possibility of investigating the glass transition in structural 
glasses from finite-size correction to mean-field theory justifying efforts in this direction. 

* * * 

We thank for useful discussions C. Donati, S. Franz, F. Sciortino and A. Rocco for a 
careful reading of the manuscript. F.R is supported by the Ministerio de Educacion y Ciencia 
in Spain (PB97-0971). 
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